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Abstract 

The classification of minimal rational surfaces and the birational links 
between them by Iskovskikh, Manin and others is a well-known subject in 
the theory of algebraic surfaces. We explain algorithms that realise links 
of type II between minimal del Pezzo surfaces, one of the major classes of 
birational links, and we describe briefly how this fits into a large project to 
implement the results of Iskovskikh's programme in Magma. 

The theory of minimal rational surfaces and their birational links works over an 
arbitrary perfect field k. Our interest here is the case /c = Q or a number field, in 
part because we can compute in these fields, but also because the implementation 
we present here imposes some conditions on the characteristic. We let k denote the 
algebraic closure of k. When a geometric object A is defined over k, A denotes its 
base change to k. All algebraic surfaces in this paper are nonsingular. A surface X 
is rational if there is an isomorphism X P| defined over k (but not necessarily 
over k). This notion is sometimes also called geometrically rational, to emphasise 
the algebraic closure in the statement. A surface X is minimal if any birational 
morphism X ^ Y defined over k is an isomorphism. 

When k = C the results of the theory are over a century old. The use of Mori 
theory and the Sarkisov programme, following Corti [C] and Iskovskikh [I], make 
our approach here very different to what it would have been 15 years ago. Yet more 
recent results, such as Hacon and McKernan's approach to the Sarkisov programme, 
may change it again in the future. 

As an introduction, we outline the results of the theory of minimal rational 
surfaces. This theory is the end result of 150 years of development: from Cayley's 
computation |Caj of the 27 lines on a cubic surface around 1850 (in effect, computing 
the divisor class group); via Castelnuovo's rationality criterion of around 1900 (which 
determines, over an algebraically closed field, whether or not the surface can be 
rationally parametrised); Segre's analysis [S^ in 1942 of the nonrationality of a cubic 
surface in terms of its Picard group, working over an arbitrary field; and Manin [M] 
and Swinnerton-Dyer's |SD] geometrical analysis of rationality in the 1960s in terms 
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of points and configurations of the 27 lines; and culminating in Iskovskikh's complete 
analysis in a series of papers during the 1980s. A full account of the modern theory 
is given in [I] ; it is also sketched in the appendix to [C] , together with details of the 
(uni)rationality of these surfaces over the ground field. The book [KSCj gives proofs 
of several parts of the theory. 

We sketch the classical and modern theories in Section [H Our main technical 
tools are explained in Section [21 and they are applied in Section [3] to compute the 
major birational maps of the theory: Geiser and Bertini involutions. In Section H] we 
show how our algorithms can be applied to analyse the Cremona group of birational 
selfmaps of a minimal cubic surface. Section [5] explains how the work described here 
fits into a broader research programme. 

We have implemented the algorithms described below in the computational al- 
gebra system Magma [BCP] : all code and examples are available at |BKR] . 



1 Minimal models of rational surfaces 

1.1 Classical results over the complex numbers 

A complete analysis of this case is in |KSC] Chapters 2 and 3. Over C, a surface is 
minimal if and only if it does not contain any curves C = with = —1. 

The minimal rational surfaces over C are: P^, the minimal rational surface scrolls 
F„ for n > 2, and Fq = P^ x P-^, which is often embedded as a quadric {xy = zt) C P^. 
Magma treats scrolls as ambient spaces using multigraded ringsj^ the scroll F„ has 
homogeneous coordinate ring k[u,v^x^y\ bi- graded by the columns of the matrix 
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and irrelevant ideal (u, v) fl (x, y), which is indicated by the separating vertical line 
in the matrix. From the point of view of Mori theory, it is crucial to consider each 
surface with a chosen structure map: for P^ it is the trivial map to a point, and 
for the scrolls F„ with n > 2 it is their natural map to P^. The quadric Fq admits 
two maps to P^, the two projections to its factors, and we must choose one of these. 
Although we often omit the map in our notation, it is illiterate not to have it (or, 
equivalently, the corresponding choice of extremal ray) in mind. 

We describe four classes of birational links, or elementary birational maps, be- 
tween these surfaces (including Fi, although it is not minimal): 

I (blowup) P^ Fl in suitable coordinates by {u,v,w) t— > {u,v,w, 1). 



-'^This is not in the 2008 export of Magma, but will be included in later versions. More generally, 
any implementation of toric geometry contains surface scrolls among its first examples. 
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II (elementary transformation) either Fj Fj_i in coordinates by {u,v,x,y) ^-^ 

{u, V, X, uy), or Fj --->• Fj+i in coordinates by (u, v, x, y) ^ {u, v, ux, y). 

III (blowdown) Fi — > in suitable coordinates by {u,v,x,y) {uy,vy,x). 

IV (change factors) the identity map on Fq, but taking the two different projection 

maps Fq on the source and target. 

The main result is that any birational map between minimal rational surfaces fac- 
torises as a composition of these elementary birational maps. 

Theorem 1.1 (Noether-Castelnuovo). Let X and Y be minimal rational surfaces 
over k = C and (f: X Y a birational map between them. Then there are 
birational links ei, . . . ,er and an automorphism ipofY such that (p = ipoSrO- ■ -oei. 

For example, if X = F = P^ and (f: {x,y,z) i— > {l/x,l/y,l/z) is the standard 
quadratic Cremona transformation then (p factorises as 

p2 Fi — ^ Fo --^ Fi P^ 

Indeed, as a composition of maps, up to a linear automorphism of P^ it is 

(x, y, z) (x, y, z, 1) (x, y, z, x) 

{x,y,z,xy) = {xz,yz,l,xy) {xz,yz,xy) = (1/y, 1/z), 

corresponding to the sequence: blowup in (0, 0, 1), elementary transform in (0, 1, 0, 1), 
elementary transform in (1, 0, 0, 1), blowdown of the negative section. 

1.2 Minimal rational surfaces over a perfect field 

We will be concerned with two classes of surfaces. First, the class of del Pezzo sur- 
faces of Picard rank 1; such surfaces are automatically minimal since any nontrivial 
(birational) morphism X ^ Y decreases the Picard rank. Second, the class of conic 
bundles over a smooth rational curve. Surfaces in this class are not necessarily min- 
imal; the surface Fi is not minimal, for instance. Nevertheless, a minimal rational 
surface over k belongs to one of these two classes. 

Minimal del Pezzo surfaces. These are surfaces X with —Kx ample and Pic(X) = 
Z. In their (pluri-)anticanonical embedding, they are in one of the families listed in 
Tabled! in this embedding, the divisor of any degree 1 linear section of X is linearly 
equivalent to -Kx- In the table, the entries X C P^ X C P^ X C P®, X C P^ 
are surfaces of degree d in P*^; these are not our main concern here, so we give only 
coarse information about their equations (although see Section [5^ for more results). 
In degree d = 5, the ideal of X C P^ is generated by the five maximal pfaffians of 
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Kx Description 

~~1 X6CP(1,1,2,3) 

2 X4 CP(1,1,1,2) 

3 X3 C P3 

4 X2,2 C 

5 X C P^ a general linear section of Grass(2, 5) C P^ 

6 X C P^ defined by 9 quadrics 

8 X C P^ defined by 20 quadrics 

9 X C P^ defined by 27 quadrics 

Table 1: Families containing the minimal del Pezzo surfaces. The case = 8 splits 
into two, according to whether X is isomorphic to a quadric X = X2 C P'^ or not. 
The case = 9 splits into two, according to whether X = P^ or not. 

a skew 5x5 matrix of linear forms — that is one interpretation of the role of the 
Grassmannian. 

We emphasise that whether a given surface X in Table [T] actually belongs to the 
class we are discussing depends on whether Pic(X) = Z; this, in turn, depends on 
both k and on the defining equations of X, and in general it is a difficult problem. 
If so then the Picard group is generated by the class of —Kx in all cases except the 
nonsingular quadric X2 C P^ (a special case for = 8), when —^Kx generates, 
and P^ (a special case for = 9), when —\Kx generates. 

Conic bundles. These are surfaces X that admit a morphism / : X — * C to a 
smooth rational curve C such that Pic(X) = /* Pic(C) +Z and a general geometric 
fibre F of / is a conic F2 C P^. The base curve C may have /c-rational points (in 
which case it is isomorphic to P^) or not (in which case it is itself isomorphic to a 
plane conic). 

If gi, . . . , ^ C are the closed points at which / is degenerate — the irregular 
values of / — then the degree of X is 

= 8 - deg(gi) deg(g„). 

When C ^ P^ the surface X can be written as a relative anticanonical model: 

X = (F = 0) in the scroll F^,,, = 

where < a < 6 and (i > are all integers, and F = F{u,v,x,y, z), in coordi- 
nates u, V, X, y, z on Fa^^ of bi-degrees given by the columns of the matrix, is bi- 
homogeneous of weight id, 2); in particular, F is a quadric in x, z. (Although only 
d > —a is required for this linear system to contain irreducible surfaces, such X 
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has a section whenever d < and so is not minimal unless X = F„; compare with 
Lecture 2 of [Rj.) In this notation, it is easy to compute 

K]^ = 8-3rf-2(a + 6), 

so J2 deg(gi) = 3d + 2(a + b). Suitable diagonal surfaces 

X = {Aax^ + Bd+2ay^ + Cd+2bZ^ = 0) C F,,^ 

with A, B, C forms in u, v of the indicated degrees give examples with arbitrarily 
complicated irregular values (the roots of ABC = where the remaining quadric is 
irreducible over k). 



1.3 Birational links 



The surfaces we consider are examples of two-dimensional Mori fibre spaces, that 
is, maps f : X —>■ S with X a surface, 5* a point or a nonsingular curve and / a 
morphism with connected fibres, —Kx relatively ample and of relative Picard rank 1. 
The map / is simply the given map for a conic bundle, and it is the trivial map 
to the point Spec/c, denoted also by {*} with k implicit, when X is a minimal del 
Pezzo surface. 

A birational link between two-dimensional Mori fibre spaces f: X ^ S and 
f : X' S' is a diagram 

— (1.1) 
/' 

S S' 
in which F is a birational map arising in one of four ways: 

Type I. These are commutative diagrams of the form 



X- 



X' 



S- 



■S' 



where F ^ is the blowup of an irreducible closed point of X. For example, 
take X = F^, S = {*} (a point), X' = Fi and S' = ¥\ 



Type II. In this case there is a surface Y and maps h: Y 
fitting into a commutative diagram: 

Y 



X and h': Y ^ X' 





X 



X' 



S- 



■S' 
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The maps h and h' are the blowups of irreducible closed points of X and X' 
respectively. 

Type III. These are inverses of links of type I, so there is a commutative diagram 

X X' 

S ^S' 

in which F is the blowup of an irreducible closed point of X'. 
Type IV. Here we have a diagram of the form 

X= X' 
S S' 




Spec k 



in which X and X' are the same surface but the link changes the Mori fibre 
space structure. For example, take X = Fq = x with X ^ S the 
projection onto the first factor and X = X' ^ S' the projection onto the 
second factor. 

The following central result is the analogue of Theorem II. 1[ 

Theorem 1.2 (Iskovskikh [I] Theorem 2.5). Let X and Y be minimal rational 
surfaces over a perfect field k and ip: X --^ Y a hirational map between them. 
Then there are hirational links ^i, . . . and an automorphism ip ofY, all defined 
over k, such that ip = o e^. o ■ ■ ■ o si. 

1.4 Geiser and Bertini involutions 

In addition to the factorisation theorem above, Iskovskikh [I] Theorem 2.6 classifies 
all the hirational links (11. ip that can occur into 41 different classes according to 
the type of X and X' . In Sections [2H51 we describe the tools central to the imple- 
mentation of seven of these classes: the Geiser and Bertini involutions on del Pezzo 
surfaces of degrees 2, 3, 4 and 5. These are all links of type II. 

We use the following notation throughout. The degree 1 hyperplane section of 
X C P*^ is denoted by A, and \nA\ denotes the linear system of all sections of 
degree n. We regard elements of \nA\ both as homogeneous polynomials on the 
whole of P'^ and as the zero loci on X they define, and we move freely between these 
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two descriptions. (When n is large, we may take a vector space complement to those 
polynomials that lie in the ideal of X — they do not define divisors on X, of course.) 
A subsystem is denoted by 7i C \nA\; this is some linear subspace of polynomials 
of degree n. If P G X is a point, then the space of polynomials of degree n whose 
zero loci on X vanish to order m at P is denoted by \nA — mP\ (or T-C{—mP) if 
restricting attention to a subsystem TC). 

Geiser involutions 

Let X C P'^ be a del Pezzo surface of degree ci = 3, 4 or 5 and P G X a closed point 
of degree d — 2. The Geiser involution ip: X --->• X with centre P is defined as 
follows: P spans a linear 11 = P'^^^, and a general hyperplane H = P'^"^ containing 
n intersects X in c? points, the sum of P and an effective 0-cycle Q of degree 2; 
ip exchanges the geometric points of Q (whether they are defined over k or not). 

This clearly defines an involution of X and it is straightforward to check that ip 
is a rational, and hence birational, map. Moreover, if we suppose that X is minimal, 
following Iskovskikh the linear system corresponding to ip is \{d — 1)A — dP\. 

Remark. The minimality condition may be surprising here; the point is that if {E C 
y) — s> (P G X) is the blowup of P G X then Riemann-Roch computes 

x{id-l)A-dE) =d+l, 

but we need to prove that this is equal to h'^(^{d — l)A — dE) , which requires the first 
cohomology to vanish. Minimality ensures this, although in practice it is often the 
case for nonminimal examples too. The same remark holds for Bertini involutions. 

Now we describe ip as a link of type II, in the style of Section [T751 We form the 
diagram 

Y 

9 

X X 

{*} W 

in which f : Y X is the blowup of P, the Mori fibre space structure X {*} 
on X is the trivial one, mapping to a point, and g: Y X is a morphism that we 
define below. 

It can be shown that there is a unique effective curve C on X that is of anti- 
canonical degree d — 2 and has multiplicity d — lai each geometric point P' G P: for 
example, ii d = 3 then C is the tangent curve Tp{X) fl X. Let F G Y he the strict 
transform of C. Then the following hold (we omit the proofs): P is a (— l)-curve or 
a union of conjugate (— l)-curves, its contraction g maps to P'^, and for some choice 
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of g — which is only defined up to an automorphism of P by what we have said so 
far — we have g: Y ^ X Gf"^ and the composite g o is equal to zp: X --■> X. 

Bertini involutions 

Let X C P*^ be a del Pezzo surface of degree c? = 2, 3, 4 or 5 and P E X a. closed 
point of degree d — 1. The Bertini involution jp : X --■> X with centre P is defined as 
follows: P spans a linear 11 = P'^"^, and a general hyperplane H = containing 
n intersects X in a nonsingular curve C of genus 1. Moreover, C has a fc-rational 
point Q (the residual point to P in X fl 11) , so (C, Q) is an elliptic curve. The map 
jp acts by —1 in the group law on this elliptic curve. 

This clearly defines an involution of X and it can be shown that jp is a rational, 
and hence birational, map. Moreover, if we suppose that X is minimal, following 
Iskovskikh the linear system corresponding to jp is \(2d — 1)A — 2dP\. As above, 
minimality is sufficient, but not necessary, for this linear system to give the right 
map. 

There is an extra detail when d = 2: to obtain a map to P^(l, 1, 1, 2), we must 
also compute sections of {{Ad — 2) A — 4dP\ and find an additional section not in the 
subspace spanned by products from \ {2d — 1)A — 2dP\. 

The Bertini involution is described as a link of type II in the same way as the 
Geiser above; the only difference is that the contracting curve F G Y, where Y is 
the blowup oi P G X, is now the unique effective curve in \{2d — 2) A — (2d — 1)E\. 

Both the Geiser and Bertini involutions are naturally defined algebraically using 
the pluri-anticanonical model of Y. This description clarifies many of the points 
we have touched on — for example, it explains why it is natural that the Bertini 
involution should act by —1 in the group law on the elliptic curve [C, Q) above — and 
it generalises to higher dimensions. See [CPRj for explanation and many examples. 

2 Imposing conditions on functions 

This is the main algorithm. Let X C P'^ be a surface and P G X a. point (/c-rational 
or not). Let Ti G \nA\ be a linear system on X that is the restriction of some linear 
system on P'^. We want to compute the linear system 7i(— mP) GTi oi divisors that 
have multiplicity at least m G N at P. 

There is one point to note. In solving these equations we compute a system of 
homogeneous polynomials on the ambient space P*^. Although we can, and do, use 
these to make a map P'^ P'^, we only work with the restriction of this to X — the 
extension to P'^ is not determined by the theory. 
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2.1 Basic case: H. = \nA\ and P E X is /c-rational 

For this discussion, we assume that X is a surface and P G X a rational point, the 
case we need for our apphcation. However, the description extends to the blowup of 
any nonsingular point on a variety of any dimension. It also applies to blowups of 
higher- dimensional centres T G X, as long as one can compute in generic coordinates 
along r. 

Normal coordinates along the blowup of P. The idea is to compute the 
blowup y of P G X as a sequence of implicit functions on a blowup patch of the 
ambient F'^. The functions defining X are denoted Fi, . . . ,Fr (that is, these form a 
basis for the ideal of X). 

Let c = d — 2 he the codimension of X C P'^, so c < r. It is easy to describe 
an affine patch ip: A'^ ^ P'^ of the blowup of P G P'^. Furthermore, we may assume 
we have coordinates ui,...,Ud on A"^ for which the exceptional divisor E d Y 
of the blowup of P G X is the Uc+i-axis (= n^^i-axis) and Ud is not a critical 
direction for the functions Ff that is, {dFj/dud)j=i^,,,^r 7^ along E. Roughly, we 
work in coordinates Ud-i,Ud on F in a neighbourhood of the generic point of the 
exceptional divisor E; precisely, this means working over the function field k{E) = 
k{ud-i) and regarding Ud as a formal power series variable. With these variables 
as coefficients, the equations defining Y are polynomials in mi, . . . ,Mc, and there is 
a basis fi,...,fc of polynomials of this ideal. We may assume, by using pairwise 
resultants if necessary, that these polynomials are each univariate in one of the 
variables: that is, fi = fi{ui) for each i = 1, . . . , c. Now the implicit function theorem 
guarantees the existence of power series fi{ud) over k{E) for which fi{(pi) = 0. These 
are computable to any given precision; Ui = ipi are the equations of Y, and we can 
use them to eliminate the variables ui, . . . ,Uc- 

Compute Ti. along the exceptional divisor. Let h be the generic element of 7i, 
a homogeneous polynomial of degree n with unknown coefficients; we may assign 
indeterminates ai, . . . , otv as these coefficients, where N is the number of monomials 
of degree d on P^. Pulling h back to the normal coordinates to on F that we 
computed above expresses it as a power series in Ud with coefficients in k{E), say 
qo+qiUd+q2Ud+- ■ ■ , computed to precision at least m™, with coefficients qi = qi{uc+i) 
rational functions in k[E) and in the unknown coefficients Oi, . . . , oat of /i. The latter 
only appear linearly in the numerator of each q^ — pulling the expression aix'^ + 
a2x'^~^y + ■ ■ ■ back to the normal coordinates evaluates the monomials x'^, x'^'^y, . . . 
at expressions in those coordinates and then gathers terms together first in powers 
of Ud and then in powers of Uc+i, the variable along E. The order of this power 
series is the degree of vanishing of Ti. along E, so sections of T-C{—mP) are those 
polynomials whose coefficients are solutions of the first m coefficients of the power 
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series, thought of as systems of hnear equations in ai, . . . , oat. 

These systems can be solved using standard computational algebra tools, and a 
basis for the solution space provides the coefficients of homogeneous polynomials of 
degree n that base 7i{—mP). 

2.2 Modifications for more general cases 

In general, we do not work with the full linear system |?t,74| but with some subsystem 
Ti C \nA\. This makes no essential difference: we use a basis of the sections of Ti in 
place of the monomials of degree n. 

The main routine must also impose conditions at points P & X oi higher degree. 
We may assume that P is irreducible over A;; it is enough then to make a finite field 
extension k G ki that splits P into /ci-geometric points and to apply the algorithm 
at just one of these. This determines /ci-linear conditions on the coefficients of the 
linear system, which in turn determine degP linear conditions over k. 

3 Del Pezzo surfaces of low degree 

To construct links of type II on del Pezzo surfaces we must compute Geiser and 
Bertini involutions given a centre P & X. We concentrate here on surfaces X of 
degree 3 or 4; we discuss surfaces of degrees 2 and 5 in Section |5] below. 

3.1 Applying the algorithm 

Suppose P G X is an irreducible closed point. The calculation is in two steps. If P 
has degree d — 2, compute a basis for the sections of \{d — l)A — dP\ to make a Geiser 
involution; if P has degree d — 1, compute a basis for the sections of \{2d—l)A — 2dP\ 
to make a Bertini involution. In either case, this basis has d+1 elements and defines 
a map tp: F'^ P*^, which (when restricted to X) is the required involution up to 
the choice of basis we made for the linear system. But without the correct basis, i/j 
is unlikely even to map X to itself. There is a unique /c-linear automorphism of P'' 
such that i)oilj is the required involution; we must find this 'missing' automorphism 

To determine we use the geometric definition of the involution. The aim is to 
find d+2 points of X that span P'' and whose images can be computed; these points, 
together with their images, are enough to determine t?. As usual, the existence of 
many fc-rational points is not expected, but intersecting X with random linear spaces 
containing P provides arbitrary numbers of closed points of low degree. 

Specifically, in the Geiser case we take general linear spaces of dimension d — 2 
through P and compute the residual intersection with X (using an ideal quotient 
to remove P from the intersection). Over the closure k, this consists of two points 
(li,Q2 G X. We can compute Qi = ip{qi) and the missing automorphism {} must 
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satisfy "^iQi) = Ri, where we set -Ri = q2- In the Bertini case we take hyperplanes 
n through P and intersect with X. If 11 is sufficiently general then X fl 11 is 
nonsingular and we can compute a Weierstrass normal form over k; general fibres 
of the elliptic involution are defined over k, so they give points of degree 2 on X 
whose components qi,q2 G X are exchanged by the Bertini involution. Again we 
can compute Qi = ip{qi) and the missing automorphism must satisfy '&{Qi) = Ri, 
where we set Ri = q2- 

It remains to solve for given d + 2 pairs of independent points Qi, Ri as above. 

3.2 Finding the missing automorphism 

Let Qi, . . . , Qd+2 and . . . , Rd+2 be two sequences of points of Pp each of which 
spans Pp and let be the /c-linear automorphism taking each Qi to Ri. Standard 
linear algebra routines compute a matrix representing -t}, but in our application 
there are two additional aspects to consider. First, our i} is defined over k so the 
corresponding matrix should have entries in k too; by itself, this is not a problem. 
Second, in practice we do not work over k but over a different extension k (Z ki G k 
for each pair of points Qi and Ri] this makes it difficult to apply the solution 
algorithm for linear equations directly. To deal with this we rephrase the algorithm 
slightly and compute a matrix representing 'd as follows. 

Let Mq be a matrix whose rows are representatives for Qi, . . . ,Qd+2 and Mji 
a similar matrix for the Rj. We seek a {d + 1) x [d + 1) matrix M for which 
MqM = Mr. It is enough to solve this up to scalar multiples of rows; so, letting 
Kr. be a matrix whose columns are a basis of ker it is enough to solve the system 
of equations 

QiMKr, =0, . . . , Qd+2MKR^_^^ = 0. (3.1) 

Each of these imposes d linear conditions on M. Since the equations are independent, 
the solution space has dimension {d + 1)^ — d{d + 2) = 1, and the entries of M (itself 
only defined up to a scalar) are then the coefficients of any nontrivial solution. This 
does not yet solve the problem: the coefficients appearing in (13.11) still lie in the 
various fields ki and it would be expensive to compute a composite field containing 
them all. 

Taking the trace of each of the equations gives a new system 

TTk,/k{QiMKR,) = 0, . . . , TTk,^^/k{Qd+2MKR^^J = (3.2) 

defined over k. The following lemma is elementary; the main point is to avoid 
/c- linear relations holding between the chosen representatives of the Qi. 

Lemma 3.1. Let Qi, . . . , Qd+2 o^nd Ri, . . . , Rd+2 be two sequences of points of P^ 
as above. For fixed representatives of the Qi and Ri, consider the systems of linear 
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equations (13. ip and (13. 2p defined onV = k^'^+^^' 



as above. Denote the solution sets 



of these two systems by WssJi and Ws3:2l respectively. 

Then W^^ji C W^jsjj and equality holds for a general choice of representative for 
each Qi. 

3.3 Degree 4 del Pezzo surfaces 

To illustrate our algorithms, we give examples to show how to construct involu- 
tions using an implementation in Magma; the implementation, together with these 
examples and their output, is available at [BKRj . 
First make a nonsingular surface 



defined over k = Q. 

> P4<x,y,z,t,s> := ProjectiveSpace(Rationals() ,4) ; 

> f := x*y - z*t + 2*x"2 + s"2; 

> g := -x~2 + y~2 - z~2 + t~2 - s"2; 

> X := Scheme (P4, [f ,g]) ; 

For a Geiser involution we need a point P G X of degree 2, which we construct by 
intersecting X with a particular line: 

> P := Intersection (X, Scheme (P4, [x,z,s])); 

> Degree (P) ; 
2 

This is all the data needed to construct the Geiser involution centred in P. The 
map is stored in Magma as a composition of simpler maps, so we use Expand (G) 
to see its defining equations of degree — 1 = 3; it is usually costly to do this step 
and is unnecessary. We also check that the involution really does map X to itself. 

> G := GeiserInvolution(X, P) ; 

> Expand (G) ; 

Mapping from: Pr j : P4 to Pr j : P4 
with equations : 

4/3*x*z'2 + 2/3*x*z*t - l/3*y*z*t - l/3*x*t'2 - l/3*x*s'2 
+ l/3*y*s~2 

-2/3*x*z~2 - y*z~2 - 7/3*x*z*t + 2/3*y*z*t + 2/3*x*t"2 

- l/3*x*s~2 - 2/3*y*s~2 
y~2*z + z*t"2 - z*s"2 

4*y~2*z - 4*z"3 - y"2*t + 4*z*t'2 - t'3 - 2*z*s'2 + t*s"2 
y"2*s + t"2*s - s"3 

> G(X) eq X; 




)4 



true 
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For a Bertini involution we need a point Q G X of degree 3, which we construct as 
the residual intersection to a 2-plane 11 containing a rational point (0, 1, 1, 0, 0) G X. 
The Support (Z) command below computes the fc-rational support of Z. 

> Pi := Scheme (P4, [x+y-z,s]); 

> Z := Intersection(X, Pi); 

> supp := Support (Z); supp; 
{(0:1:1:0:0)} 

> Degree (Z) ; 
4 

> R := Cluster (Representative (supp) ) ; 

> Q := Scheme ( P4, ColonIdeal(Ideal(Z) , Ideal(R) ) ); 

> B := BertiniInvolution(X,Q) ; 

Computing the map B takes about half a minute, giving a map defined by polyno- 
mials of degree 2d — 1 = 7\ the equations of the map are large, and so here we only 
show the initial terms. 

> B; 

Mapping from: Pr j : P4 to Pr j : P4 

Composition of Mapping from: Pr j : P4 to Pr j : P4 

with equations : 

y~2*z''5 - 301231/288*y~2*z~4*t - 102767/144*x*z~5*t 

- 11755/32*y*z-5*t + 101059/72*z-6*t 

- 4791269/4608*y"2*z-3*t-2 - 2205985/2304*x*z~4*t~2 - . . . 

> B(X) eq X; 
true 

These maps are indeed involutions and it is not necessary to check this explicitly — 
although one can check, for instance, that G o G is the identity map by using the 
interpolation routines explained below. 

3.4 Examples on a cubic surface 

First make the surface X: (x^ + 2y^ + 3^;^ + = 0) C defined over = Q. 

> P3<x,y,z,t> := ProjectiveSpace(Rationals() , 3); 

> X := Scheme(P3, x"3 + 2*y"3 + 3*z"3 + 4*t"3) ; 

To make a Geiser involution, we need to choose a rational point for its centre. We do 
not have code for finding rational points — in general, this is an unsolved problem — 
but in this case we can see some obvious choices: (1,— 1,-1, l)eX, for instance. 

> p := X ! [1,-1,-1,1] ; 

> G := GeiserInvolution(X,p) ; 
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For a Bertini involution, we need a point of degree 2. The line L = {y + z + t = 
x — z + t = 0) meets X in the rational point (3, 1, 1, —2), and the residual intersection 
to this is an irreducible point of degree 2. 

> L := Scheme ( X, [y+z+t, x-z+t] ); 

> Z := [ Y : Y in IrreducibleComponents (L) | Degree (Y) eq 2 ] [1] ; 

> B := BertiniInvolution(X,Z) ; 

(The Bertini calculation takes around 4 seconds; by comparison, the Geiser calcu- 
lation is instant.) The map B is defined by fairly large polynomials of degree 5. 
We can construct another birational selfmap by composing these. 

> h := Expand (B * G) ; h; 
Mapping from: Prj : P3 to Prj : P3 
with equations : 

54518131/19784704*x-4*y'6 + 59844679/7419264*x~3*y~7 + 
272174051/29677056*x-2*y-8 + 7725367/1391112*x*y-9 + 
9154681/5564448*y-10 + 82450383/9892352*x~4*y'5*z + ... 

The equations of h are large polynomials of degree 10 that run over several pages. 

4 The group of birational selfmaps of a cubic sur- 
face 

It is well known, [M] Theorem 38.1 for instance, that Geiser and Bertini involutions 
(together with the subgroup of linear automorphisms) generate the group Bir(X) of 
birational selfmaps of any minimal cubic surface X = X^ C P^. This also follows 
from Iskovskikh's classification of links [T] Theorem 2.6: the only elementary links 
from cubic surfaces are birational selfmaps. (This is not the case in degrees 4 and 5 
where factorisations of birational selfmaps into elementary links may pass through 
other surfaces.) 

The proof of Theorem 11.21 works by induction on the degree of the given self- 
map (p. The idea is to find a basepoint of degree 1 or 2 that has high multiplicity 
in curves belonging to the linear system defining if. Given such a point, a so-called 
maximal centre, the proof precomposes (p by the Geiser or Bertini involution it de- 
termines; this decreases the degree and the induction continues. We implement 
this algorithmic step, demonstrating it here by factorising the map h computed in 
Section [331 

First find the base locus of h. This is done naively, setting the defining equations 
of h to be zero. The result could be strictly bigger than the base locus, but since 
we check the multiplicity of base points later this does not matter. 

> base_h := Scheme (X.DefiningEquations (h) ) ; 
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> Diiiiension(base_h) ; 


> Degree (base_h) ; 
205 

We need to identify irreducible components of this base locus, and this could be a 
problem if there really were 205 base points. But of course the calculation above 
has found a non-reduced scheme, and since we only need to know the base points 
set theoretically, we can reduce it before further analysis. Unfortunately, this seems 
to be difficult; it takes about 4 minutes. 

> base_red := ReducedSubscheme (base_h) ; 

> Degree (base_red) ; 
3 

Since the map h is not linear, there must be a maximal centre in this base locus. Un- 
twisting by Bertini involutions is likely to reduce the degree of h more dramatically 
than by Geiser involutions, so we look for maximal centres of degree 2 first. 

> Q := [ Y : Y in IrreducibleComponents (base_red) | Degree (Y) eq 2 ] [1] ; 

Of course, Q is exactly the enforced base point Z from section lH^ above. but without 
knowing that we need to check that it is a maximal centre. 

> assert is_maximal_centre(X, Q, h) ; 

This takes about a minute: it runs the main algorithm in high degree to check that 
Q has multiplicity strictly greater than 10 in the linear system defining h. We could 
omit this step since checking that the degree of h is reduced after untwisting is 
sufficient. 

We must now untwist h by the Bertini involution centred in Q. (Notice the order 
of composition in MAGMA: this is hi = ho ei.) 

> epsl := BertiniInvolution(X, Q) ; 

> hi := epsl * h; 

There is a practical computational point here: if we expanded out the equations 
of hi, they would have degree 10 x 5, the product of the degrees of the equations 
of El and h. But untwisting is meant to reduce the degree. Working modulo the 
equation of X, the equations of hi have a large common factor that can be cancelled, 
but it is not clear how to do this calculation. 

Instead we use interpolation to find the the correct equations. Bertini involutions 
reduce the degree by 4(multQ(/i) — deg(/i)), so even without knowing multQ(/i), the 
possible degrees for the resulting equations are limited. Given a target degree, the 
interpolation evaluates hi at many points of X and uses this collection of points and 
their images to impose linear conditions on the coefficients of the desired equations. 
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As usual, we do not have a supply of rational points of X to work with (there may 
be none), but we can intersect X with random rational lines to get points of degree 3 
and use these. If the target degree was too low, the solution space will only include 
multiples of the defining equation of X. The first time there are other solutions, 
these will be the coefficients of the map. The function we use below returns a 
boolean value, which is false unless there is a unique additional solution. If the 
boolean value is true, the function also returns the equations of hi. (The assert 
bool statement causes a crash unless the boolean is true.) 

> bool,eqns := interpolate (X, hi ,2) ; 

> assert bool; 

We rebuild hi with the lower-degree equations. 

> hi := map< P3 -> P3 I eqns >; hi; 
Mapping from: Prj : P3 to Prj : P3 
with equations : 

x*y + y~2 + 3/2*x*z + 3/2*z~2 + 2*x*t - 2*t"2 

l/2*x~2 + l/2*x*y + 3/2*y*z - 3/2*z'2 + 2*y*t + 2*t~2 

l/2*x~2 - y"2 + l/2*x*z + y*z + 2*z*t + 2*t"2 
-l/2*x"2 + y"2 + 3/2*z"2 + l/2*x*t + y*t + 3/2*z*t 

Repeat the process with hi. 

> base_hl := Scheme(X, Def iningEquations (hi) ) ; 

> Dimension(base_hl) ; 


> basel_red := ReducedSubscheme (base_hl) ; 

> Degree (basel_red) ; 
1 

> r := Representative (Support (basel_red) ) ; r; 
(1 : -1 : -1 : 1) 

So the base locus is a single rational point (1,-1,-1,1) G X. There is no choice 
but to untwist by the Geiser involution. 

> eps2 := GeiserInvolution(X, r) ; 

> h2 := eps2 * hi; 

> bool, eqns := interpolate (X,h2 , 1) ; 

> assert bool; 

> h2 := niap< P3 -> P3 | eqns >; h2; 
Mapping from: Prj : P3 to Prj : P3 
with equations : x, y, z, t 

The resulting map /12 gives a linear automorphism of X so the factorisation is com- 
plete. In this case, it is easy to check that the identity is the only automorphism 
of X, so the group Bir(X) is generated by all Geiser and Bertini involutions — what 
this group is, therefore, has become an arithmetic question involving the low degree 
points of X and the relations between Geiser and Bertini involutions. 
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5 Rational surfaces and MAGMA 



We describe briefly the broader context of a possible complete implementation of 
minimal rational surfaces in Magma that builds on the algorithms here. It is realis- 
tic to expect: to compute the factorisation of rational maps between two-dimensional 
Mori fibre spaces, the so-called Sarkisov programme; to address rationality questions; 
and to analyse elliptic fibrations. Beck and Schicho have algorithms, implemented 
in Magma, that compute the Picard group of a del Pezzo surface and, if necessary, 
carry out the minimal model programme in certain circumstances. 

5.1 The Sarkisov programme for rational surfaces 

We were able to carry out the Sarkisov programme for selfmaps of a cubic sur- 
face in Section H] because in that case other surfaces are not involved. The same 
is true for del Pezzo surfaces of degrees 1 and 2, but in other cases we need the 
full classification of links in order to proceed. An implementation of the full Sark- 
isov programme for factorising birational maps between minimal rational surfaces, 
following Iskovskikh — so realising Theorem 11.21 explicitly — is feasible, given the fol- 
lowing additional components. 

The involutions on del Pezzo surfaces of degrees 2 and 5 are essentially the 
same as those described here. In degree 5 the calculations are currently too slow to 
be practical. In degree 2 we must work with weighted projective space; this only 
complicates the calculation slightly. The main additional difficulty is in recognising 
the surface as an image of the multi-linear system that determines the involution. For 
example, if P = (1, 0, 0, 0) G X C P(l, 1, 1, 2) with tangent plane Tp{X) = {y = 0), 
in coordinates x, y, z, t, then 

\3A-AP\ = {y^yh,yz',y^x) 

and |6A— 8P| is spanned by quadratic expressions in these together with yH and y^zt. 
These determine a map X ---> P^(l^,2^) whose image Y is isomorphic to X — but 
we would need to compute | — Ky\ and | — 2Ky\ to make that identification (and 
then study the geometrical description of the Bertini involution as before to make 
the right choice of automorphism). 

Computing automorphisms of del Pezzo surfaces is straightforward. Since del 
Pezzo surfaces are embedded by the anticanonical class, any automorphism X —* X 
extends to an automorphism of the ambient projective space. If this space is P'^ 
then the only automorphisms are projective linear maps; when X = X^ C wP^ = 
P^(l, 1, 1, 2) we must allow quasi-linear maps as follows: 

xi^ fi{xi,X2,X'i),. . . , X3 /3(xi,X2,X3), y ^ ay + g{xi, X2, x^) 

where xi,X2,xs,y are the homogeneous coordinates on w¥^, the /, are linear forms, 
(? is a quadratic form and a & k does not have a square root in k. 
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General conic bundles can be embedded in scrolls, even if the base rational curve 
does not have a rational point. With such a description, one can compute links by 
explicit equations as we have for del Pezzo surfaces. Cremona and van Hoeij [CvHj 
give an algorithm for Tsen's theorem (over an extension of k, as necessary), which 
is one of the essential tools for this. 

5.2 Rationality questions for minimal rational surfaces 

While we assume that our surfaces X are rational over k, the question of whether 
they are also rational over k is well studied. Swinnerton-Dyer |SDj answers the 
question for a cubic surface in terms of a rational point on X and conditions on the 
configuration of the 27 lines. In degree 5, Enriques proved that X automatically 
has a rational point and is rational over k; see |Sh-B] . In degrees > 6, X is rational 
over k if and only if it has a rational point, and either case can happen. In degrees 
8 and 9, de Graaf-Harrison-Pilnikova-Schicho [dG et al] determine whether a del 
Pezzo surface (not necessarily minimal, and for k a number field) is rational and 
they compute a parametrisation over k in that case; this is implemented in Magma. 

5.3 Elliptic fibrations on minimal rational surfaces 

Elliptic fibrations birational to minimal rational surfaces are classified in some cases: 
Dolgachev [D] classifies elliptic fibrations birational to (the construction of these 
fibrations dates back to Halphen [H] ) ; Cheltsov |Ch] and Brown-Ryder |BR] analyse 
elliptic fibrations birational to minimal cubic surfaces; and Cheltsov also analyses 
the case of del Pezzo surfaces of degrees 1 and 2. A birational map from X to an 
elliptic fibration can be regarded as a limiting case of a birational map to another 
Mori fibre space, so the methods of construction and exclusion for the two problems 
are very similar. The computational aspects of the elliptic fibration problem are 
fully analysed in the degree 3 case in |BR] . and a full Magma implementation is 
given; this is not yet done in the other cases. 
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